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Abstract 


Light  scattering  measurements  were  carried  out  to  determine  the  aerosol 
properties  of  the  atmosphere.  First  the  expression  for  the  aerosol  attenuation 
coefficient  is  derived,  based  on  the  field  geometry  in  conjunction  with  Rayleigh 
and  aerosol  scattering  considerations.  Then  the  results  derived  from  the 
measurements  are  discussed.  The  paper  concludes  with  an  atlas  of  105  profiles 
for  altitudes  to  about  35  km  (the  data  does  not  exclude  the  presence  of  aerosols 
with  low  number  density  between  35  to  70  km).  These  profiles  consist  of  aerosol 
attenuation  and  extinction  coefficients  as  a  function  of  altitude.  Since  the  co¬ 
efficients  are  proportional  to  aerosol  number  density,  the  profiles  yield  informa¬ 
tion  concerning  aerosol  stratification.  A  plot  for  the  computed  mean  of  the  105 
vertical  profiles  is  included. 
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Symbols 


2 

A  =  area  (normal  to  beam  axis)  of  scattering  volume  (cm  ) 
B  =  base  line  distance  between  projector  and  collector  (cm) 
D  =  diameter  of  projector  mirror  (cm) 

o 

E  =  irradiance  (watts /cm  ) 

=  instrumentation  response  (volts) 

Fd  =  total  energy  at  detector  (watts) 

F  =  total  energy  of  source  (watts) 
h  =  altitude  (km) 

J (4  )  =  intensity  of  angular  scattered  energy  (watts/ster) 
s 

L  =  length  of  scattering  volume  (cm) 

.  i 

P  (<t>  )  ~  normalized  aerosol  phase  function  (sr  ) 

P  8 

P  J  )  =  normalized  Rayleigh  phase  function  (sr  *) 
r  s 

R  =  distance  between  source  and  scattering  volume  (cm) 

=  total  transmission  between  source  and  scatter  volume 

Tp  =  aerosol  transmission  for  both  slant  paths 

Tr  =  Rayleigh  transmission  for  both  slant  paths 

T  *  total  transmission  between  source  and  scatter  volume 
s 

Z  =  zenith  angle 


vii 


A  =  aerosol  attenuation  coeffic  lent  (km  ) 

P  .} 

4  =  Ravleigh  attenuation  coefficient  (kin  ) 

r  - 1  - 1. 

3(ij>  )  =  totai  volume  scattering  function  (cm  ster  ) 

0  =  beam  divergence 

X  =  wavelength)  (m) 

t  =  aerosol  optical  thickness 
P 

b  =  angle  between  beam  and  collector  field 

^  =  elevation  angle  of  detector  mirror 

6  =  elevation  angle  of  projector  mirror 

P 

<*>  =  scattering  angle 

‘  ,  solid  angle  subtended  by  collector  with  apes  a.  element  ot  scattering  volume 
(sterachan) 


t 
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An  Atlas  o •  Aerosol  Attenuation  and  Extinction  Profiles 
for  the  Troposphere  and  Stratosphere 


I.  INTRODUCTION 

In  nearly  all  problems  concerned  with  the  interaction  of  light  with  the  atmos¬ 
phere,  the  aerosol  attenuation  coe  ficient  in  the  troposphere  and  stratosphere 
emerges  as  an  indispensible  and  1  itle  known  parameter.  But  it  is  an  important 
parameter  because  it  describes  quantitatively  how  the  aerosol  content  of  the  atmos¬ 
phere  effects  the  passage  of  light.  Also,  because  of  its  proportionality  to  the 
aerosol  number  density,  this  parameter  provides  information  concerning  aerosol 
stratification.  There  has  been  much  emphasis  on  aircraft,  balloons,  and  rockets 
in  the  pursuit  of  this  research  disc  ipline.  However,  ground-based  remote  sensing 
methods  have  great  merit.  First,  they  have  the  important  capability  of  performing 
true  in  situ  measurements.  Then  here  is  the  added  advantage  of  being  able  to 
acquire  large  quantities  of  data  at  relatively  small  er  st. 

In  this  report,  the  searchlight  method  of  optical  probing  used  in  New  Mexico 
is  summarized.  Although  measurement  capability  at  0.  55  micron  wavelength 
exceeds  70  km,  the  altitudes  of  interest  will  be  discussed  up  to  35  km  because 
above  this  altitude  the  measurements  show  that  aerosol  attenuation  for  most  pur¬ 
poses  can  be  neglected.  A  description  of  the  methods  used  was  published  recently 
(Elterman,  1966).  At  this  time,  t!  e  primary  objective  is  to  present  an  atlas  of 
profiles. 


(Received  for  publication  28  November  1966) 


Table  1.  Computer  Output  (partial  tabulation) 

Measurement  on  13  April  1364  at  00:18 


2.  ENERGY  TRANSFER  EQUATIONS 

Figure  1  presents  the  searchlight  scene  geometry.  The  projecting  mirror  is 
parabolic  with  a  diameter  of  91.  .  cm,  and  its  elevation  angle  is  fixed.  The  collec¬ 
ting  mirror  is  identical,  ant’  sens  up  or  down  the  beam  with  elevation  angle  <f> 

The  beam  and  field  of  view  are  t  ligned  precisely  so  as  to  be  co-planar  for  all 
elevation  angles.  Approximately  30  minuter  is  required  for  a  continuous  scan  of 
the  beam  over  the  altitude  range  '.  8  to  35.  3  km  (0  to  57°  elevation).  The  altitude 
of  the  scattering  volume  is  established  by  the  into:  section  of  the  optical  axis  of  the 
beam  and  the  collector  field  of  view .  The  altitude  resolution  averages  about  1  km. 
A  partial  list  of  the  geometrical  relationships  is  shown  in  Table  1. 
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Since  the  divergences  of  the  beam  and  collectin'*  field  are  small,  the  opposing 
extremities  of  the  scattering  volume  are  nearly  equal,  both  angularly  and  in  length. 
From  *he  scene  geometry  (see  Figure  1),  the  equations  for  the  energy  transfer 
are: 


(1)  Source  to  Scattering  Volume 


F  •  T  *  EA 
s  s 


(1) 


where: 


Fg  =  total  energy  of  source  (watts) 

Tg  »  transmission  ^molecular  +  aerosol)  source  to  scatter  volume 

A  «  area  (normal  to  beam  axis)  at  scattering  volume,  receiving  incident 
2 

radiation  (cm  ) 

o 

E  *  irradiance  incident  on  area  A  (watts/cm  ) 


(2)  At  the  Scattering  Volume 


J  (6  )  «EALp  (<!>) 

8  o 


(2) 


where:' 


<t>s  *  scattering  angle 

J($g)  *  intensity  of  angularly  scattered  light  (watts/steradian) 

(<f> s )  a  (molecular  +  aerosol)  angular  scattering  function  (cm  *ster  ) 

L  =  length  of  scattering  volume  (cm) 


(3)  Scattering  Volume  to  Detector 

Fd*J(V  TdW 


(3) 


where: 

F^  *  total  energy  at  detector  (watts) 

T  ,  =  transmission  (molecular  +  aerosol)  scatter  volume  to  detector 
d 

w  =  solid  angle  subtended  by  collector  with  apex  at  element  of  scattering 
volume  (steridiars) 

Combining  Eqs.  (1),  (2)  and  (3),  the  energy  at  the  detector  is 


rd"F.  TsTdL“lJ|*s> 


M) 
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From  Figure  1  it  follows  that  the  product  Lw  is  constant;  that  is. 


L  «  R  0/ sin 

^B 

(a) 

Sin  *B  “  f 

sin  </>p 

(6) 

*D2 

w  *  - 

(7) 

4R 


and  combining  Eqs.  (5),  (6)  and  (7), 


Lw 


*  D 2  0 
4B  sin  <t> 

P 


(8) 


Here  *  75°  and  C,  is  the  constant. 

The  instrumentation  response  to  the  scattered  energy  is  linear  so  that 
Erp  =  Cg  Fd<  Also,  measurements  show  that  the  variation  of  source  intensity  for 
each  scan  is  small  so  that  F  can  be  assigned  the  constant  C».  If  the  constants 
are  combined,  C  =  C,,  C^,  and  Eq.  (4)  becomes 

Erp  =  CTsV<V  (9) 


Rearranging  the  transmissions  and  expressing  the  volume  scattering  function 
by  its  components, 

Ts  Ta  .  Tr  TP  (10) 


and 


W*.>+*PpP 

s 

'  where; 

Tg  *  total  transmission  source  to  scatter  volume 
Td  =  total  transmission  detector  to  scatter  volume 
Tr  =  Rayleigh  transmission  for  both  slan‘,  paths 
Tp  =  Aerosol  transmission  for  both  slant  paths 


(11) 


c 


* 


K 


Rayleigh  attenuation  coefficient  (cm  *) 
Aerosol  attenuation  coefficient  (cm  *) 


-U 


P  to„)  =  normalized  Rayleigh  phase  function  (steradian  ) 

4  8 

P  to  )  =  normalized  aerosol  phase  function  (steradian  *) 
p  s 

Equation  (9)  takes  the  form 


Erp  =  C  Tr  Tp  [fir  •  Pr  tog)  +  /3p  •  Pp  tos)J 


(12) 


where: 

Erp  =  instrumentation  response  (volts) 

C  =  proportionality  constant  (volt  cm  steradian) 


3.  AEROSOL  SCATTERING  EQUATION  AND  ITS  SOLUTION 


Equation  (12)  has  the  limitation  that  the  constant  C  is  unknown.  The  procedure 
for  evaluating  C  is  demonstrated  by  referring  to  Figure  2  which  represents  a  set 
of  measurements  obtained  8  May  1964  at  03:05.  The  ordinate  is  designated  as  the 
relative  response  since  these  values  are  normalized  to  unity  at  35.  3  km.  Also 
shown  is  the  calculated  Rayleigh  response  for  identical  conditions.  As  in  other 
measurements,  the  curves  nearly  converge  over  the  last  few  kilometers. 
Accordingly,  at  about  35  km,  aerosol  scattering  is  assumed  to  be  sufficiently 
small  so  that  it  can  be  neglected.  Then  Eq.  (12)  for  these  conditions  becomes 

Erp(35)  =  C  Tr(35)  T  (35)  0r<35)  Pr(132°)  (13) 

Here  132°  corresponds  to  the  scattering  angle  for  the  35.3  km  altitude.  Com¬ 
bining  Eqs.  (12)  and  (13),  the  altitude  variation  of  the  aerosol  attenuation  coefficient 
is  expressed  by 

E  (h)  T  (35)  T  (35)  P  (132°) 

a  (h)=_£E - £ - E - s  os) _ i _ 

P  Erp(35)  Tp00  Tp(h>  Pp»s) 


-  Br(h) 


P„  J 
r  s 

P  to  J 
P  s 


(14) 


RELATIVE 


0  5  10  15  20  25  30  35 

ALTITUDE  (km) 


Figure  2.  Instrumentation  Response  of  Single  Scan  rnri  Corresponding 
Rayleigh  \  alues;  Wavelength  is  0.  55  Microns 
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As  mentioned,  the  terms  E _ (h)  and  E  (35)  are  obtained  from  the  measure- 

rp  rp 

ments.  All  the  Rayleigh  terms  are  calculated  by  using  the  Atmospheric  Attenuation 
Model  tabulations  (Elterman,  1964)  for  0.55m  wavelength.  The  Rayleigh  phase 

o 

Assuming  for  the  present  that  the  aerosol 
phase  function  is  established,  then  Eq.  (14)  has  one  unknown  because  T^(h)  can  be 
expressed  in  terms  of  (h) 


function  is  Pr  (</>g) 


0.75  (J  +  cos*  *g). 


Tp(h) 


exp 


ii 

-Z 


2.76 


1 3p  (h)  A  h  sec  Z(h) 


(15) 


where:' 


/Ip  =  mean  value  of  P  for  the  atmospheric  layer  (km  *) 

Ah  *  layer  thickness  (km) 

sec  Z(h)  =  sum  of  the  secants  of  the  zenith  angles  for  the  beam  and  collector 
field  slant  paths 


Since  /3  exists  in  the  transmission  term  as  an  exponent,  the  solution  of 

Eq.  (14)  is  not  analytic.  The  method  used  is  an  iterative-convergent  procedure. 

The  iteration  is  started  by  using  the  Atmospheric  Attenuation  Model  (Elterman, 

1964)  for  0.  55  u  to  calculate  an  initial  series  of  aerosol  transmission  values, 

Tp  (h),  applicable  to  the  field  geometry.  Having  done  this,  the  Model  values  of 

Tp(h)  have  no  further  function.  The  iteration  proceeds  as  shown  in  Figure  3  until 

convergence  is  attained.  Equation  (1  i)  then  yields  a  final  set  of  P  (h). 

P 

As  to  the  choice  of  the  aerosol  phase  function,  P  (A  ),  the  literature  is  ex- 

P  8 

tensive  with  both  computational  and  experimental  studies.  Here  must  be  con¬ 
sidered  the  work  of  Waldram  (1945),  Reeger  and  Siendentopf  (1946),  Bullrich  and 
Moeller  (1947  and  1964),  Foitzik  and  Zschaeck  (1952),  Deirmendjian  (1964), 
Fraser  (I960),  Barteneva  (1960)  and  othei  s  who  made  important  contributions  to 
this  area  of  research.  The  criteria  finally  considered  for  the  choice  of  the  phase 
function  were  (a)  its  derivaiion  should  be  from  field  measurements  with  a 
meteorological  range  of  at  least  30  km,  and  (b)  the  phase  function  should  permit 
Eq.  (14)  to  converge  within  a  reasonable  number  of  iterations.  The  results  of 
Reeger  and  Siedentopf  (and  there  are  others)  readily  meet  these  criteria.  The 
aerosol  phase  function  was  derived  from  their  measurements  (sec  Figure  4). 


NORMALIZED  phase  function 
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Figure  3.  Iteration  Procedure.  Tp  is  aerosol  transmission. 
ftp  is  aerosol  attenuation  coefficient 


Figure  i.  Normalized  Phase  Functions 
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4.  I\;>TIU  MENTATION 


The  methods  follow  those  used  previously  (Elterman  1962,  and  Elterman  et  al 
1964),  but  with  considerable  modification.  The  beam  intensity  is  modulated  by  a 
shutter  mechanism  fronting  the  searchlight.  Figure  5.  The  modulation,  20 cycles/ 
second,  represents  a  compromise  between  a  desired  high  modulation  frequency 
and  the  mechanical  spin  limitation  of  the  shutters.  An  auxiliary  detector,  mounted 
on  the  beam  modulator,  monitors  the  source  intensity  and  generates  a  signal 
which  is  synchronous  with  the  modulation.  The  synchronous  signal  is  transmitted 
to  a  s.te  30.2  km  distant,  where  the  optical  collector,  synchronous  detector  and 
amplifier  are  located.  The  collector  mirror  and  the  photomultiplier  assembly 
comprise  she  sensing  device.  The  detector-filter  response  peaks  at  0.55  microns 
and  has  a  wavelength  bandwidth  at  the  half-amplitude  response  approximating  550 
Angstroms. 


5.  SKY  B'VCMiROl.ND 

The  measurements  were  conducted  during  moonless  nights  so  that  the  sky 
background  would  be  minimized.  Although  modulation  of  the  beam  contributes 
greatly  to  isolating  the  desired  scattering  from  the  sky  background  and  other 
spurious  sources,  it  does  not  completely  effect  isolation.  The  sources  of  the 
background  level  are:  (a)  the  star  field  and  airglow  (direct  and  scattered), 

(b)  secondary  and  higher  order  scattering  from  the  beam,  and  (c)  instrumenta¬ 
tion  noise. 

Probably  the  major  source  is  direct  and  scattered  light  from  the  star  field 
and  the  night  airglow.  Beam  modulation  with  synchronous  rectification  is 
effective  m  suppressing  these  sources,  but  not  entirely  so.  The  incidence  of 
modulated  and  unmodulated  eneigy  on  the  photomultiplier  results  in  a  shot  noise 
spectrum,  part  of  which  is  passed  by  the  bandwidth  of  the  amplifier  (0.02  cycles/ 
second).  Secondary  and  higher  order  scattering  from  the  beam  enters  the  field 
of  the  collector  mirror  for  the  most  part  at  low  altitudes  (regions  of  relatively 
high  atmospheric  molecular  density),  but  is  considered  small  due  to  the  30.2  km 
separation  of  the  sites.  Finally  the  photomultiplier  dark  current  and  other 
instrumentation  noise  contribute  to  the  signal  background.  The  signal-to-notse 
ratio  at  35  km  exceeds  20,  and  increases  rapidly  at  lower  altitudes. 

b.  OBSKRY  \ T10NS  YM)  RESIT, TS 

The  data  obtained  is  in  the  form  of  a  continuous  trace  on  the  chart  recorder. 
The  elevation  angle,  is  recorded  auto;  atically  on  the  chart  along  with  the 
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Figure  5.  Modulating  Shutters  with  Drive  Motors  Below.  Detector 
assembly  for  synchronous  rectification  mounted  near  the  top 


response  to  the  scattered  light.  Subsequently,  the  response  is  digitized  at  1° 
intervals  as  in  Table  1.  The  results  based  on  sampling  of  a  year's  measurements 
show  that  it  is  convenient  to  categorize  the  aerosol  distributions  derived  from  the 
data,  into  "moderate-structured",  "mediur'-structured",  and  "full-structured" 
profiles. 

Following  through  on  the  measurements  discussed  in  Figure  2,  the  response 
values  were  used  in  Eq.  (14),  and  the  resuiting  aerosol  attenuaiion  coefficients 
were  plotted  automatically  from  the  punch  card  compute  r  output  as  shown  in 
Figure  6.  The  Rayleigh  coefficients  and  the  extinction  coefficients  (Rayleigh  + 
aerosol)  are  included.  The  character  of  the  aerosol  curve  places  it  in  the 
category  of  a  "medium-structured"  profile.  Since  the  aerosol  coefficients  pri¬ 
marily  are  proportional  to  the  aerosol  number  density,  examples  of  these  cate¬ 
gories  will  be  discussed.  The  diminution  of  the  aerosol  density  values  over  the 
lowest  few  kilometers  occurs  frequently.  A  region  of  aerosol  c  oncentration  is 
evident  with  a  maximum  at  18  krn.  For  the  region  2.  8  km  down  to  sea  level. 


8  MAY  1964 
03505 


••• 


RAYLEIGH  (COMPUTED) 

-f-  +■  +-  AEROSOL  (FROM  MEASUREMENTS  WITH  SMOOTHING ) 

. EXTINCTION  (RAYLEIGH  B  AEROSOL) 

- EXTRAPOLATED  (IN  ACCORDANCE  WITH  ELTERMAN,  1964) 


l _ I - 1 - i— - 1 - 1 - 1— 

o  5  10  15  20  25  30 

ALTITUDE  (km) 

Figure  6.  Aerosol  and  Related  Parameters.  The  aerosol  profile 
represents  the  "medium-structured"  category 
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extrapolation  was  carried  out  to  conform  with  the  0.  55  n  tabulations  of  the  Atmos¬ 
pheric  Attenuation  Model  (Elterman,  19G4). 

The  profiles  in  Figure  7  are  examples  from  the  "moderate -structure" 
category.  Here  the  aerosol  number  densities  between  2.76  and  about  14  km  show 
relatively  little  change.  The  identity  of  the  20  km  layer  is  fairly  evident. 

The  staggered  scale  presentation  in  Figure  8  shows  a  series  of  profiles 
taken  over  a  six-hcur  period.  These  profiles  are  categorized  as  "full -structured". 
This  body  of  data  was  obtained  consecutively  except  for  periods  of  interim  activity 
pertaining  to  operational  requirements.  The  low  altitude  concentration  shows  a 
maximum  at  5  km  hi  profiles  A,  B,  and  C.  During  the  night,  the  maximum 
becomes  less  defined  and  in  profile  E  acquires  a  simple  logarithmic  distribution 
(2.  8  to  6.  5  km).  Next  a  tropopause  concentration  with  a  maximum  at  about  9  km 
is  evident,  but  this  region  increasingly  loses  its  identity  in  profiles  D  and  E. 

There  follows  a  region  of  aerosol  concentration  above  the  tropopause  (12  to  23  km) 
which  persists  through  the  night,  and  is  identified  with  the  20  km  aerosol  layer 
determined  by  the  direct  measurements  of  Junge,  Chagnon  and  Manson  (1961). 
Finally,  a  26  km  maximum  snows  surprising  temporal  integrity.  Above  28  km, 
the  decreasing  values  indicate  a  quasi -Rayleigh  atmosphere. 

The  usefulness  of  the  0p(h)  profile  for  obtaining  quantitative  information  can 
be  readily  demonstrated.  In  Figure  8  profile  C,  the  aerosol  concentration 
usually  designated  as  the  20  km  layer  can  be  examined  with  considerable  detail. 

It  is  well-defined,  asymmetrical  and  extends  in  altitude  from  12.7  to  23.2  km 

-4 

(Table  1).  The  median  of  the  layer  occurs  at  18. 2  km  where  8  -  9.  33  x  10 

- 1  P 

km  and  the  turbidity  factory,  8/8  =0.  83.  The  maximum  aerosol  density, 

P  r  _3  . 2 

however,  occurs  at  15.  6  km  where  8  s  1.66  x  10  km  and  the  turbidity 

1^ 

factor  is  0.  99.  The  aerosol  optical  thickness  of  the  layer  is  determined  by 
23.2 

-rp=^  Jp(h)Ah  (16) 

12.7 


where: 


■1. 


8  (h)  is  the  mean  aerosol  attenuation  coefficient  (km  )  for  each  altitude 
P 

increment 

Ah  is  tne  altitude  increment  (km)  as  designated  in  Table  1, 


This  leads  to  the  value  rp  =  0.01  whicn  is  modest  but  important  in  applications 
requiring  large  zenith  angles.  The  slant  path  transmissions  through  the  layer  due 
to  Rayleigh  and  aerosol  attenuation  were  obtained  with  Eq.  (15)  and  are  shown  in 
Figure  9. 
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Figure  7.  Aerosol  Profiles  Representing  the  "Moderate-Structured 
Category 


AEROSOL  ATTENUATION  COEFFICIENTS  (km) 
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Figure  8.  Aerosol  Profiles  Representing  the  "Full-Structured1 
Category  Showing  Variations  Over  a  Period  of  About  Six  Hours 
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7.  KRKOR  CONSIDERATIONS 

Here,  it  will  be  convenient  to  designate  a  set  of  measured  values,  E  (h),  as 
the  input  profile  and  the  computed  aerosol  attenuation  coeffi  nents,  (3  (h),  as  the 

r 

output  profile. 

Examination  of  Eq.  (141  shows  that  the  term,  Erp  (35),  can  bear  significantly 
on  the  accuracy  of  the  results.  It  is  necessary,  therefore,  to  consider  the  effect 
of  noise  and  background  on  this  term.  If  the  background  level  causes  Erp(35) 
to  increase,  then  some  negative  values  in  the  output  profiles  may  result.  Also, 
the  measured  and  Rayleigh  curves  (as  shown  in  Figure  2)  are  nearly  converged 
between  30  and  35  Ion.  The  presence  of  noise  then  may  cause  the  measured  value 
to  be  less  than  the  corresponding  Rayleigh  value,  which  is  an  unacceptable  condi¬ 
tion.  Accordingly,  the  procedure  was  adopted  where  the  input  profile  is  smoothed 
for  the  last  few  kilometers.  The  computer  then  examines  the  derived  £Sp(h)  values 
after  the  first  iteration  for  the  presence  of  negative  values.  Routinely,  the 

instrumentation  is  set  so  that  the  measured  value  of  E  (35)  is  near  unity.  If  for 

rp 

h<35  km  any  negative  values  of  0p  are  present,  Erp(35)  is  reduced  by  0.001  re¬ 
peatedly  until  cnly  positive  values  of  0  (h)  exist.  The  iteration  then  continues  with 
no  further  adjustments  until  convergence  is  obtained.  This  procedure  corrects 
Ef  (35)  for  noise  and  background. 

The  angular  pointing  error,  when  scanning  the  beam  upward  or  downward, 
approximates  ±  0.  1  deg.  ee.  This  is  equivalent  to  an  error  of  0.  05  to  0.  09  km, 
depending  on  the  altitudes  of  measurement. 

The  February  1964  data  above  22.6  km  may  be  questionable  due  to  uncertain 
operating  procedure. 

It  is  reasonable  to  estimate  that  the  other  functional  aspects  of  the  measure¬ 
ments  can  result  in  about  5  or  6  percent  error  in  the  Erp(h)  values.  If  systematic 
(always  negative  or  positive),  then  this  percentage  applies  to  the  entire  input 
profile,  and  the  resulting  output  profile  will  be  unaffected  because  we  are  dealing 
with  ratios,  Erp(h)/Erp(35).  However,  a  part  of  the  input  profile  subjected  to 
a  5  percent  change  results  in  considerable  error  due  to  the  summation  procedure, 
Eq.  (15).  This  is  demonstrated  by  choosing  a  severe  example  of  error  simulation;' 
that  is,  increasing  the  Erp(h)  values  over  a  substantial  portion  of  the  input  profile 
(approximately  27  km  of  Table  1).  If  these  are  increased  by  5  percent,  the  output 
profile  values  of  /3p  change  28  percent  at  2.  76  kin  and  40  percent  at  29.7  km.  A 
corresponding  6  percent  increase  results  in  changes  of  (3p  from  33  to  48  percent. 

In  instances  where  convergence  cannot  be  attained,  the  data  is  abandoned. 
Usually  failure  to  converge  is  due  to  high  values  of  Erp(h)  for  altitudes  under  5  km, 
caused  by  heavy  dust  or  thin  cloua 
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Five  major  factors  determined  the  number  of  profiles  acquired  during  the 
period 

(1)  The  operational  time  that  could  be  allowed  was  limited. 

(2)  Operations  were  conducted  only  during  moonless  nights  in  order  to  keep 
noise  background  at  a  minimum. 

(3)  Operations  were  cancelled  when  meteorological  conditions  were  unfavor¬ 
able.  This  included  high  winds  as  well  as  cloudy  skies.  July  and  August  season¬ 
ally  are  months  when  both  dust  storms  and  cloudy  skies  occur  daily  in  New 
Mexico. 

(4)  Equipment  difficulties  accounted  for  absence  of  data  in  January  1964  and 
at  other  times. 

(5)  Failure  of  site  facilities,  especially  telephone  lines,  resulted  ir.  con¬ 
siderable  loss  of  time.  Synchronous  rectification  and  amplification  depended  on 
transmittal  of  the  synchronizing  signal  by  telephone. 

During  the  period  12  December  1963  through  12  December  1964,  a  total  of 
105  useful  profiles  were  acquired,  listed  as  follows:. 


Month 

Profiles 

Month 

Profiles 

Dec  1963 

22 

July  1964 

3 

Jan  1964 

Aug  1964 

-- 

Feb  1964 

12 

Sept  1964 

2 

Mar  1964 

6 

Oct  1964 

6 

April  1964 

28 

Nov  1964 

4 

May  1964 

9 

Dec  1964 

5 

June  1964 

8 

The  atlas  of  profiles  that  follow  suggests  that  characteristic  stratifications 
exist.  Although  these  are  designated  as  layers,  they  do  not  appear  to  be  discrete. 
Rather,  the  profiles  indicate  an  aerosol  continuum  having  greater  or  lesser  con¬ 
centrations.  Although  results  are  presented  for  altitudes  up  to  about  35  km,  the 
data  acquired  does  not  exclude  the  presence  of  aerosol  with  low  number  density 
between  35  and  70  km.  A  computed  mean  of  105  profiles  comprising  the  atlas  is 
shown  in  Figure  10. 

As  mentioned,  the  measurements  were  conducted  at  0.  55  microns.  If  for 
exploratory  calculations,  aerosol  attenuation  coefficients  are  desired  for  other 
wavelengths,  one  can  consider  using  the  results  in  conjunction  with  an  empirical 
relationship  such  as  ~  X  . 

The  salient  purpose  of  this  paper  is  to  make  a  substantial  number  of  profiles 
available,  possibly  for  use  in  various  disciplines.  At  this  time  it  is  left  to  the 
reader  to  cull,  tabulate,  compare,  and  contemplate  this  body  of  material.  We 
plan  to  continue  evaluation  and  further  study  of  these  results. 
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Atlas  of  Profiles 

This  section  contains  a  series  of  105  useful  profiles,  for  altitudes  to  about 
35  km,  acquired  during  the  period  12  December  1963  through  12  December  1964. 
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U.  ABSTRACT 


Light  scattering  measurements  were  carried  out  to  determine  the  aerosol 
properties  of  the  atmosphere.  First  the  expression  for  the  aerosol  attenuation 
coefficient  is  derived,  based  on  the  field  geometry  in  conjunction  with  Rayleigh 
and  aerosol  scittcring  considerations.  Then  the  results  derived  from  the  measure¬ 
ments  are  discussed.  The  paper  concludes  wit.,  an  atlas  of  105  profiles  for 
altitudes  to  about  35  km  (the  data  does  not  exclude  the  presence  of  aerosols  with 
low  number  density  between  35  to  70  km).  These  prof*1'®  consist  of  aerosol 
attenuation  and  extinction  coefficients  as  a  function  of  a  •  iae.  Since  the 
coefficients  are  proportional  to  aerosol  number  density,  me  profiles  yield  informa¬ 
tion  concerning  aerosol  stratification.  A  plot  for  the  computed  mean  of  the  105 
verticai  profiles  is  included.  ,  -~ 
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